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Abstract 

Background: Grapevine {Vitis vinifera L) is an economically innportant fruit crop. Quality-deternnining grape connponents 
such as sugars, acids, flavors, anthocyanins, tannins, etc., accumulate in the different grape berry development stages. 
Thus, correlating the proteomic profiles with the biochemical and physiological changes occurring in grape is of 
paramount importance to advance in our understanding of berry development and ripening processes. 

Results: We report the developmental analysis of Vitis vinifera cv. Muscat Hamburg berries at the protein level from fruit 
set to full ripening. An iTRAQ-based bottom-up proteomic approach followed by tandem mass spectrometry led to the 
identification and quantitation of 41 1 and 630 proteins in the green and ripening phases, respectively. Two key points in 
development relating to changes in protein level were detected: end of the first growth period (7 mm-to-15 mm) and 
onset of ripening (15 mm-to-VlOO, VlOO-to-1 10). A functional analysis was performed using the Blast2G0 software based 
on the enrichment of GO terms during berry growth. 

Conclusions: The study of the proteome contributes to decipher the biological processes and metabolic pathways 
involved in the development and quality traits of fruit and its derived products. These findings lie mainly in metabolism 
and storage of sugars and malate, energy-related pathways such as respiration, photosynthesis and fermentation, and the 
synthesis of polyphenolics as major secondary metabolites in grape berry. In addition, some key steps in carbohydrate 
and malate metabolism have been identified in this study, i.e., PFP-PFK or SuSy-INV switches among others, which may 
influence the final sugar and acid balance in ripe fruit. In conclusion, some proteins not reported to date have been 
detected to be deregulated in specific tissues and developmental stages, leading to formulate new hypotheses on the 
metabolic processes underlying grape berry development. These results open up new lines to decipher the processes 
controlling grape berry development and ripening. 
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Background (FAO, 2007 [1]). In the food industry, grape berries are 

Grapevine {Vitis vinifera L.) is the most important fruit widely commercialized raw, as a fresh and dried fruit, or 

crop in the world, and viticulture and enology play an are processed mainly as grape juice or wine. Extracts from 

important role in the economy of many developed and grape berry seeds, skin and also grape plant leaves are used 

emerging countries. In 2007, the Food and Agriculture as commodities in nutraceutical and cosmetic industries. 

Organization (FAO) of the United Nations estimated that The grape berry is a non climacteric fruit that exhibits a 

the gross world production of grape berries was over 67 double sigmoid growth pattern [2]. The first growth phase 

million tons per year, with land occupancy of 7,272,583 ha after fruit set is characterized by not only rapid cell div- 
ision, which increases the number of cells, but also by an 
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characterized by important biochemical and physiological 
changes such as softening, coloring and engustment of 
berry. As grape berries develop, they change in size and 
composition. These transformations range from small, 
firm and acidic with little sugar, desirable flavors or aroma, 
to becoming larger, softened, sweet, highly flavored, less 
acidic and highly colored fruit. The biochemical changes 
underlying grape berry ripening have been recently 
reviewed [3]. The development of these traits is subject to 
genetic, environmental and viticultural conditions, which 
determine the final quality of the berry and its derived 
products. During the first growth period, chlorophyll is 
the main pigment present in fruit and cells are rich in or- 
ganic acids; the most prevalent compounds are tartaric 
and malic acids, which accumulate mainly in skin and 
flesh. Compounds vastly increase starting at veraison, and 
the major ones are glucose and fructose [4], in addition to 
phenolic and aromatic compounds [5], while malate con- 
centration lowers [6]. Flavor development in grapes is 
partly due to the acid/sugar balance [7], which is particu- 
larly important in table grapes. Muscat Hamburg is a clas- 
sical cultivar of black table grape that is grown in many 
parts of Europe, and is greatly appreciated for its pleasant 
Muscat flavour [8,9]. 

Correct grape ripening is fundamental for both commer- 
cial fruit use and wine quality. For this reason, grape berry 
development has been investigated by monitoring the ex- 
pression of quality trait-related genes in a classic targeted 
approach [4,10-13], and also by global mRNA expression 
profiling using cDNA or oligonucleotide microarrays 
[14-19] since genomic resources for Vitis vinifera and re- 
lated species have been generated on a large scale. 

In recent years, proteomics -based technologies have been 
successfully applied to grapevine to analyze stress responses 
[20-23] and berry development [24-28]. These studies have 
improved the understanding of changes in grape berries at 
the proteome level, but are limited by the intrinsic inaccur- 
acy of single staining-based quantitative two-dimensional 
gel electrophoresis (2-DE) [29]. The application of more ro- 
bust quantitative proteomic approaches, such as difference 
gel electrophoresis (DIGE) or MS -based stable isotopic la- 
beling, should lead to a better coverage and accuracy of 
protein changes. DIGE has been applied to study grape 
berry development. It has led to a better coverage and ac- 
curacy of the protein changes reported from a gel-based ap- 
proach [30,31]. Liicker et al [32] studied grape berry 
ripening initiation by the quantitative MS-based iTRAQ 
approach, and reported a significantly larger number of 
total protein identification and protein changes than single- 
stain 2-DE-based studies. These authors noticed the diffi- 
culty of protein identification as no finished genome se- 
quence data for grapevine were available. Thus, they built a 
protein database by combining in-house generated [33] and 
publicly available grapevine nucleotide sequences. To date. 



this limitation has been greatly reduced thanks to grapevine 
genome sequencing initiatives [34,35], which have released 
the translated ORFs found in the genome to public protein 
databases. Nevertheless, the extent of sequence description 
and annotation of grapevine proteins, including genome- 
derived ones, is still far from complete. The exocarp of 
grape berries has been analyzed throughout ripening by the 
iTRAQ technology to offer better proteome coverage [36] 
and by making full use of recent grapevine genome sequen- 
cing initiatives [34,35]. 

Stable isotopic labeling followed by a reverse phase high- 
performance liquid chromatography coupled with tandem 
mass spectrometry (HPLC-MS/MS) analysis is a relatively 
recent quantitative technique, and was first described by 
Ross et al [37]. Isobaric tags for relative and absolute pro- 
tein quantitation (iTRAQ) are currently being used suc- 
cessfully to characterize and quantify changes in protein 
levels in complex biological samples [38]. The use of 
iTRAQ has become a consolidated technique in quantita- 
tive proteomics since large fold changes of protein expres- 
sion within broad dynamic ranges of protein abundance 
can be measured quite accurately [39]. The iTRAQ 
method allows the multiplexed identification and quantita- 
tion of proteins in four different samples (4-plex) and has 
been recently scaled up to measure protein changes in up 
to eight different samples (8-plex) [40]. 

The tissue-specific analysis of grape berry is critical if in- 
formation on its individual contribution to final fruit qual- 
ity is sought. Quality traits relating to the organoleptic 
properties of berry and wine are linked to specific berry 
tissue: skin composition plays an important role in deter- 
mining color, aroma, astringency and bitterness, while the 
pulp contributes to the sugar/acid balance, crispiness, 
juiciness and alcoholic potential. The seed also contributes 
to the astringency and bitterness of wines. 

Grimplet et al [41], used a large-scale transcriptomic 
analysis to study the differential expression of genes 
among main tissues, seed, flesh and skin, and reported 
that 60% of these genes exhibited a significant differential 
expression in at least one of the three major tissue types. 
Several proteomic analyses at the tissue-specific level have 
been carried out [24,26,31,32,36,42]. Nevertheless on the 
whole, a gap still remains either because one only tissue is 
analyzed or berry development is not comprehensively 
covered. 

In this study, the issue of protein expression changes 
is addressed in the pericarp/mesocarp tissues in grape 
berry and all the developmental phases are covered. For 
the purpose of gaining a better understanding of the 
protein composition in grape berry and its changes dur- 
ing development, seven different stages from fruit set to 
fully ripened were analyzed. The aim is to improve 
coverage and to complement the proteomic analysis 
done in previous DIGE experiments [31] with the high- 
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throughput proteomic analysis presented herein using 
an alternative approach, namely the top-down quantita- 
tive proteomics iTRAQ technique. The use of annota- 
tion and analysis tools [43], along with recent literature 
findings, has also led to gain insights into the biological 
significance of the proteomic data in the search for use- 
ful protein biomarkers to track the final quality traits 
associated with grape berry development. 

Results and discussion 

Grape berry development was analyzed by the iTRAQ 
technique. One iTRAQ experiment covered four time 
points of the first growth period, which spanned the stages 
from fruit set to pre-veraison (FS-to-15 mm). Another 
iTRAQ experiment covered three time points of the 



second growth period from onset to full ripening (VIOO- 
to-140g/l), and a pre-veraison fourth time point (15 mm), 
shared with the first experiment. Figure 1 describes the ex- 
perimental design and workflow for both experiments. 

Protein identification and quantitation 

In the green stages, 542 proteins were identified, which is 
similar to the number of proteins identified with at least 
two peptides in an iTRAQ experiment of the grape berry 
exocarp [36] (Additional file 1). Of these, 524 proteins 
were quantified and, after applying a p<0.05 filter, the list 
shortened to 411 (Additional file 2). All the 524 quantified 
sequences were assigned a description and 508 were anno- 
tated with an average number of five annotations per se- 
quence and an average annotation level depth of 5.8. 
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Figure 1 Workflow of the iTRAQ experiments for both the first growth period and grape berry ripening. The first growth period, 
analyzed in the first iTRAQ experiment, covers four developmental stages: fruit set (FS), 4 mm, 7 mm and 15 mm labeled with the 114, 115, 116 
and 1 17 iTRAQ tags, respectively. The ripening period covers three developmental stages, 100% veraison (V-100), 1 10 g/1 and 140 g/1, and a fourth 
point, the 15 mm stage, which is common with the previous experiment, and was co-analyzed in a second iTRAQ experiment. For each 
experiment, the extracted proteins were trypsin-digested and peptides were labeled with iTRAQ tags. After labeling, peptides were pooled and 
fractionated by strong cation exchange (SCX) chromatography. Each fraction was resolved by reverse phase chromatography and was analyzed 
directly by mass spectrometry. Searches and quantitation analysis were carried out using the ProteinPilot vl.O software (Applied Biosystems). 
Finally, the identified proteins in the set were annotated and functionally analyzed using the Blast2G0 v2.4.0 application [43] based on gene 
ontology (GO) terms. 
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In the ripe stages, 1117 proteins were identified, which 
is almost twice the number of those identified in the ex- 
periment for green stages (Additional file 3). Of these, 
1034 proteins were quantified and, after applying a 
p<0.05 filter, the list shortened to 630 (Additional file 4). 
Of the 1034 sequences, 1010 were assigned a description 
and 977 were annotated with an average number of five 
annotations per sequence and an average annotation 
level depth of 5.5. 

As a summary of the annotation results, pie charts 
for the biological function and cellular component 
gene ontology (GO) terms are shown for each experi- 
ment in Additional file 5, respectively. The resulting 
annotation files are provided as Additional file 6 
(green) and Additional file 7 (ripe). 

Proteome changes during development 

To better characterize the most active stages as regards 
proteome changes, the relative abundance of the dereg- 
ulated proteins was computed between contiguous stages 
(reporter ion ratios 115/114, 116/115 and 117/116 per ex- 
periment) and these values were used to rank the selected 
protein lists. For each stage transition, two subsets of pro- 
teins were established: those whose fold change was above 
1.5 and those with a fold change under -1.5. Figure 2 pro- 
vides a global picture of the proteome changes between 
each pair of contiguous sampled time points of the grape 
berry formation phase. Figure 2A illustrates the number 
of up- and down-regulated proteins passing the filter, 
whereas Figure 2B depicts the average and the highest 
fold change of that set of proteins. Of the number of 
deregulated proteins, three classes of development stages 
can be formed: stages with very few changes (transitions 
from 4 mm-to-7 mm and from 110 g/l-to-140 g/1); stages 
with many changes (transitions from 7 mm-to-15 mm and 
from 15 mm-to-VlOO); stages with an intermediate num- 
ber of changes (transitions FS-to-4 mm and VIOO- 
to-110 g/1). If the amplitude of changes was considered, 
the average fold change among the paired stages was not 
so different. However, there were few very major changes 
at certain transitions that revealed potential marker pro- 
teins for the developmental stages. The results obtained 
herein agree with those reported in a previous proteomic 
experiment, which analyzed the same time points using 
the DIGE platform [31]. The additional time point in- 
cluded in the present study, the transition from VlOO-to- 
15 mm, resulted in the development stage where the lar- 
gest number of proteins underwent changes and where 
the largest fold-changes occurred. These results suggest 
that the most important changes at the protein level took 
place at the end of green period and at the onset of ripen- 
ing, as previous data from genomic [44], oligo/microarray 
transcriptomic [14,15] and proteomic [27,31,32] studies 
have also indicated. 



Functional analysis 

The present study contributes to significantly increase the 
coverage of the differential proteome of the whole devel- 
opment of grape berries and confidently allows the detec- 
tion of putative marker proteins based on a previous 
proteomic experiment, which revealed that the protein 
changes observed were developmentally associated with 
the analyzed biological individuals [31]. The good agree- 
ment between most of the profiles of the proteins identi- 
fied in common herein and in Martinez-Esteso et al [31], 
serves as validation the present study using the non gel ac- 
curate quantitative proteomic iTRAQ approach. 

Functional annotation analysis 

A Fischers enrichment analysis [45] was carried out by 
comparing the frequencies of the annotation terms in the 
up- and down-regulated protein subsets with those of the 
whole list of quantified proteins as reference list. As a re- 
sult, a series of GO terms, which were statistically overrep- 
resented in the subset, was obtained (p<0.005). Figure 3 
shows the result of the Fischers enrichment analysis for 
the molecular function (F term) (Figure 3A), biological 
processes (P term) (Figure 3B) and the cellular compo- 
nents (C term) (Figure 3C) -selected GO terms whose 
level in most cases was 4 or deeper. Additional file 8 pro- 
vides a list of the proteins annotated with each enriched 
GO term. 

The proteins annotated with the GO-terms enriched 
from the up-regulated list of proteins in FS-to-4 mm were 
mainly different isoforms of the histones, tubulins and pro- 
teins involved in cytoskeleton organization and biogenesis. 
Several ribosomal proteins, which were annotated with the 
terms Intracellular nonmembrane-bound organelle-GO: 
0019685' (C term) and cellular component organization 
and biogenesis- GO:0016043' (P term), appeared to be up- 
regulated in FS-to-4 mm. In contrast, these terms were 
down-regulated from 7 mm onwards, showing those as 
leading processes during the early stages of development, 
when fruit cells were very active in division and expansion, 
and new proteins were being synthesized. The term 'hy- 
drolase activity, acting on glycosyl bonds-GO:0016798' 
(F term), which appeared to be significantly up-regulated, 
related to different glycosyl hydrolases and an acid vacuolar 
invertase, which started to accumulate from the 4 mm 
stage and was possibly involved in the carbohydrate catab- 
olism of sucrose unloaded in fruit by acting as a sink organ 
in the vine. The proteins involved in 'flavonoid biosynthetic 
processes-GO:0009813' and aromatic amino acid synthesis- 
GO:0006725' (P term) appeared to be down-regulated at 
this point of development, including the phenylalanine 
metabolic process-GO:0006558' GO term (P term), a pre- 
cursor of all the phenylpropanoid-derived compounds. This 
finding reflects an accumulation of flavonoids in early 
fruit development. The 4 mm-to-7 mm transition was 
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characterized by a few terms enriched in the down- and 
up-subsets of proteins, which reflect that most of the pro- 
cesses involved in development occurred earlier than or 
after this time point. Only the term photosynthesis light 
harvesting-GO:0009765' was up-regulated and some terms 
annotating several isoforms of histones were down- 
regulated Cnucleosome-GO:0000786' (C term) protein- 
DNA complex assembly-GO:0065004' (P term) and 'DNA 
binding-GO:0003677' (F term)), which displayed greater ac- 
cumulation earlier in grape berry development. In the next 
transition (7 mm-to-15 mm), only the down-regulated sub- 
set of proteins contained GO-enriched terms. The most 
remarkable processes, these being the 'macromolecule 
biosynthetic process-GO:0009059' (P term) and cellular 
component organization and biogenesis- GO:0016043' 
(C term), were linked with proteins like elongation factors, 
tubulins, histones and ribosomal proteins, which agrees 
with the state of cells at this point of development. These 
results tie in with the fact that fruit entered a lag phase, 
with little or no division or expansion, and most processes 
were switched off before triggering the veraison. Another 
important down-regulated term was the phenylpropanoid 
metabolism (G 0:0009699), involving flavonoids and amino 
acid precursors, which included the enzymes of the general 



phenylpropanoid pathway. The processes down-regulated 
at 7 mm-to-15 mm occurred mainly in the cytosol, unlike 
those down-regulated from 15 mm-to-VlOO, and which 
also occurred mainly in the chloroplast. Several proteins 
belonging to the photosynthetic machinery ('Photosystem 
I-GO:0009522' (C term), 'Photosystem II-GO:0009523' 
(C term)) and those involved in photosynthesis, both 
light (GO:0019684) and dark reactions (GO:0019685) 
(phosphoribulokinase, rubisco), were strongly down- 
regulated before veraison was triggered. Some related 
processes occurring in the chloroplast such as starch 
metabolism (4-alpha-glucanotransferase), redox state 
(thioredoxin, peroxiredoxin) and protein folding (chap- 
eronin 60), are linked with the terms enriched in the 
down-regulated subset, thus also indicating a decrease 
of the photosynthetic-dependent processes. In parallel, 
a significant number of ribosomal proteins was down- 
regulated (the term structural molecule activity-GO: 
0005198' (F term)). The term 'integral to membrane-GO: 
0016021' (C term) included an aquaporin and plasma 
membrane H^-ATPase, indicating a decrease in water- 
associated specific transporters prior to veraison. Likewise, 
cell wall-associated processes (endo-xyloglucan tranferase, 
pectinesterase) appeared to be deregulated, indicating that 
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(See figure on previous page.) 

Figure 3 Enriched GO terms (Molecular function (F term) (A), Biological process (P term) (B) and Cellular component (C term) (C)) for 
the sequences annotated in the up- or down-regulated subsets of proteins in each grape berry developmental stage analyzed. Bar 

diagrams indicate tine number of sequences tliat were up- and down-regulated (X-axis), annotated witli eacli enriclied GO term (Y-axis) in tine 
consecutive paired stages. GO terms descriptions are detailed in Additional file 8. Color code for the consecutive pair stages: FS-to-4 mm 
transitions (dark green bars), 4 mm-to-7 mm transitions (medium intensity green), 15 mm-to-7 mm (light green), VlOO-to-15 mm transitions 
(yellow), 110 g/l-to-15 mm transitions (light purple), 140 g/l-to-1 10 g/l (dark purple). 



these profiles play a relevant role at veraison as this stage 
was characterized by fruit softening. This is a complex 
regulated process as inferred by the terms deregulated 
at this transition, 15 mm-to-VlOO, in which some spe- 
cific isoforms related to cell wall biogenesis (the terms 
cytoplasmic membrane-bound vesicle-GO:0016023'(C term), 
'extracellular region-GO:0005576' (C term), endo-xyloglucan 
transglycosylase, expansin, polygalacturonase, pectin methyl- 
esterase, invertase pectin methylesterase inhibitor family) 
were strongly up-regulated, while other isoforms were 
down-regulated at this point. Interestingly, several 
stress-related proteins annotated with the term cop- 
per ion binding-GO:0005507' (F term) were down- 
regulated (polyphenol oxidase, diphenol oxidase, super- 
oxide dismutase). The terms which up-regulated at veraison 
were cytoplasmic membrane-bound vesicle-GO:0016023', 
extracellular region-GO:0005576' and, apart from the 
proteins indicated above, they involved ripening-related 
proteins (grip22), plasma membrane proteins and some 
proteases (subtilisin, cysteine proteinase inhibitor). Never- 
theless, the most relevant terms during this period 
were cell wall metabolic process-GO:0010382' (P term), 
chitin catabolic process-GO:0006032' (P term), invol- 
ving defense proteins (PR-4, chitinase IV), carbohydrate- 
related processes, which included the proteins involved 
in sucrose metabolism (sucrose synthase, sucrose- 
phosphate synthase) and, behind the term lyase activity- 
GO:0016829' (F term), hexose and pyruvate metabolism 
(enolase, aldolase, pyruvate decarboxylase) was found. 
Altogether, this reflects the intense activity around 
carbohydrate metabolism at the onset of ripening, pre- 
sumably in the biosynthetic direction, to convert existing 
pools of malate into sugars for their storage since gly- 
colysis is supposed to be inhibited after the onset of 
ripening [6]. 

In the VlOO-to-110 transition, the down-regulated 
process predominated the up-regulated ones. Among the 
up-regulated ones, the terms chitin catabolic process' and 
polysaccharide metabolic process-GO:0005976' (P term), 
including defense proteins, were kept, and new terms ap- 
pear such as prefolding complex-GO:0016272' (C term) 
(needed for the stabilization of the newly synthetized pro- 
teins) and cysteine-type endopeptidase activity-GO: 
0004197' (F term) (cysteine proteinase, the cysteine 
protease component of the protease-inhibitor complex 
and cathepsin b-like cysteine protease). Among the 



down-regulated proteins, the most remarkable feature 
was the diminished abundance of proteins synthesis, as 
interpreted by the enriched terms Intracellular non- 
membrane-bound organelle', (C term) cellular component 
organization and biogenesis' (P term), which included an 
important representation of different ribosomal proteins, 
and the term 'translation-GO:0006412' (P term) (t-RNA 
synthetases). Apart from these, other terms of the cellular 
component and molecular function categories associated 
with the ribosomal function were enriched in the down- 
regulated subset (macromolecule biosynthetic process' 
(P term), structural constituent of ribosome'-GO:0003735 
(C term)). Another remarkable feature was that the proteins 
involved in the photosynthetic dark reaction were strongly 
down-regulated, which reflects that grape berries lose their 
photosynthetic capability as they ripen. Finally in the 110- 
to-140 transition, only the term nutrient reservoir activity- 
GO:0045735' (F term) was enriched in the up-regulated 
subset of proteins. This term included three seed storage 
proteins that are considerably accumulated. 

Protein profile analysis 

As shown in Additional file 9, the behaviour of the time 
course of the ftinctional protein groups, which lie beneath 
important quality traits of the pericarp/mesocarp during 
fruit development, was observed in detail. The ftinctional 
protein groups of 'transporters', polyphenols', photosyn- 
thesis, respiration and fermentation' and carbohydrate and 
malate metabolism' are discussed herein. The remaining 
ftinctional protein groups, including nitrogen and amino 
acid (Additional flle 10) 'terpenoid metabolism', 'signaling 
and hormone', 'stress', 'protein synthesis', 'protein degra- 
dation', 'protein processing', 'cell division and growth, 
biogenesis' 'defense proteins' and 'other proteins of inte- 
rest', are discussed in Additional flle 11. 

7. Transporters 

Here we will focus on transport phenomena involved in 
tonoplast and plasma membrane energization for its special 
relevance in sugar transport and accumulation. The vacuole 
is a very important organelle in grape berry cells because of 
its dominant cytoplasm occupancy [46], its predominant 
role as a storage compartment of organic acids, sugars, sec- 
ondary metabolites and amino acids, and because of its 
central cell flanction in cytosolic pH regulation. To date, 
knowledge of the transporter proteins in the tonoplast is 
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surprisingly scarce. Thus, the high acidity of grapes and the 
changes that occur during ripening still mean that import- 
ant questions are raised. V-ATPase and V-PPase are two 
tonoplast primary proton pumps that are widely distributed 
in plants cells. Here, the protein levels for both pumps were 
quantified in the mesocarp during ripening. Several V- 
ATPase subunits were found up-regulated during develop- 
ment starting at the onset of ripening (Figure 4, !')♦ This 
agrees with the profiles of subunits C and F reported in 



previous studies [31]. The protein level of V-PPase proved 
constant from the onset of ripening to 100% veraison, but 
then decreased until the end of ripening (Figure 4, '2'), which 
agrees with reported transcript and protein profiles [47]. On 
the basis of V-PPase and V- ATPase proton pumping activ- 
ities, a predominant role of V-PPase in the energization of 
the grape tonoplast during development has been proposed 
despite its poorer stoichometry efficiency as compared to V- 
ATPase [47] while the role of V- ATPase became relevant 
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highlighted. The protein levels of the regulated enzymes are shown in coloured squares, indicating the change of expression (log2 ratio) for each 
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throughout the ripening phase [48]. Thus, the protein pro- 
files found herein for both the H"^ pumps were in accord- 
ance with their reported activities [48]. Although V-PPase 
was not detected during the first growth period, the profiles 
found showed more abundance at veraison than at ripening; 
such decrease combined with the emergency of a putative 
cytosolic sPPase (Figure 4, '3') along ripening, points to pyro- 
phosphate (PPi) as an energy source to pump protons across 
the tonoplast during green development, but not during 
ripening. 

Furthermore, the induction of ripening is known to be as- 
sociated with the shift of phloem unloading from the sym- 
plastic to the apoplastic pathway [28]. The latter involves 
the participation of several plasma-membrane sugar trans- 
porters which are induced at veraison [3]. Sugar 
translocation is a secondary transport driven by a proton 
gradient established by plasma membrane H^-ATPase 
(PM-ATPase). Here the abundance of PM-ATPase 
(Additional file 9A) peaked at the onset of ripening and 
became down-regulated in the following stages, which 
is consistent with previous findings [49]. This profile is 
consistent with the sugar concentrations in the phloem 
and berry apoplast; that is, high at veraison and early 
ripening, but lower as the berry ripens [50], and it re- 
sponds to the PM energization need. Thus, PM-ATPase 
can play a role relevant in triggering sugar accumula- 
tion in berries. Unlike PM, the energy demand at the 
tonoplast level continuously increased for use in both 
sugar translocation to the vacuole and the maintenance 
of the high sugar gradient concentration between the 
cytosol and the vacuole. Indeed the induction of proton- 
sugars antiporters supported vacuolar sugar accumulation 
[3]. The specific tonoplast antiporter carrier for hexoses, 
VvHT6, was up-regulated during ripening (Additional file 
9A) in accordance with reported transcript levels [14,17]. 
In addition, vacuolar H"^ pumps must play a key role in 
generating a proton motive force to accumulate sugars. As 
discussed above, both V-PPase and V-ATPase can play this 
role, but considering their respective profiles and the 
progressive accumulation of soluble PPase, the importance 
of the ATP-driven pump increases as the berry ripens. 

As note above, PPi may play key functions during berry 
development although its metabolism is still a largely un- 
known process. PPi is released as a by-product during sev- 
eral anabolic reactions, which are highly active during 
young development, such as the synthesis of DNA, RNA, 
proteins, carbohydrates, etc., and it can be simply hydro- 
lyzed by soluble pyrophosphatase (sPPase) to pull up bio- 
synthetic reactions or used in other processes. According 
to the abundance profiles of the PPi-utilizing cytosolic en- 
zymes (Figure 4), in the developing grape berry it might 
act as an energy donor per se [51] for tonoplast energiza- 
tion through the V-PPase pump, but also as a phosphoryl 
group donor in sucrose degradation via sucrose synthase 



(SuSy) (Figure 4, '4'), and UDP-glucose pyrophosphorylase 
(UDP-GP) (Figure 4, '5') and in glycolysis through 
pyrophosphate-dependent phosphofructokinase (PFP) 
(Figure 4, '6') (reviewed Plaxton, [52]). Besides the pro- 
files for these enzymes, it is interesting to note that the 
profile for sPPase correlates positively to that of phos- 
phofructokinase 1 (PFK-1) (Figure 4,7') and V-ATPase, 
but negatively to that of PPi-utilizing enzymes such as 
PFP and V-PPase, therefore, sPPase is the candidate to 
control the cytosolic PPi pools in the ripening berry de- 
velopment phase. 

2. Polyphenols 

From a wine quality perspective, the major phenolic 
compounds in grape berries are non flavonoid 
hydroxycinnamic acids (HCyAs), and flavonoid anthocy- 
anins (AGs) and proanthocyanidins (PACs) [53]. The 
synthesis of phenolics share the phenylpropanoid path- 
way, from phenylalanine (Phe) to p-coumaric acid, from 
which several branches derive, leading in grapes to lig- 
nins, hydroxycinnamic acids, stilbenoids and flavonoids 
[3] (see Figure 5). The obtained results clearly show en- 
hanced phenylpropanoid synthesis from very early devel- 
opment as the two first enzymes of the pathway, 
phenylalanine ammonia lyase (PAL, Figure 5 '1') and 
cinnamate 4-hydroxylase (C4H, Figure 5 '2'), were at 
least 4-fold up-regulated toward fruit set in relation to 
15 mm. Likewise, five enzymes of the shikimate pathway 
(3-deoxy-D-arabino heptulosonate 7-phosphate synthase 
(DHAPS, Figure 5 '3'), dehydroquinate dehydratase- 
shikimate dehydrogenase (DD-SDH, Figure 5 '4'), 
5-enol-pyruvilshikimate-phosphate synthase (EPSPS, Fig- 
ure 5 '5') and chorismate synthase (CHOSy, Figure 5 
'6')), leading to aromatic amino acids, had the same pro- 
file as the phenylpropanoid enzymes during green devel- 
opment (see Figure 5). Flavonoids constituted a 
significant portion of phenolics and included AGs, flavo- 
nols, flavan-3-ols and PACs or condensed tannins (poly- 
mers of flavan-3-ols and flavan-3,4-diols), and were the 
most abundant class of soluble polyphenolics in grape 
berries. Common enzymes for the flavonoid pathway 
(chalcone synthase (CHS, Figure 5 7'), chalcone isomer- 
ase (CHI, Figure 5 '8'), flavonol-3-hydroxylase (F3H, Fig- 
ure 5 '9'), dihydroflavonol reductase (DFR, Figure 5 10'), 
anthocyanidin synthase/leucoanthocyanidin dioxigenase 
(ANS/LDOX, Figure 5 11') and the branching point to- 
ward PACs biosynthesis (anthocyanidin reductase (ANR, 
Figure 5 12')) displayed the same abundance pattern as 
the shikimate and phenylpropanoid enzymes. Although 
no enzyme of stilbene or hydroxycinnamic acid synthesis 
branches was detected, a resveratrol/hydroxycinnamic 
acid O-glucosyltransferase (RHCA-GT) with an abun- 
dance profile like that above was identified. It shared a 
94% homology with a bi-functional glucosyltransferase 
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from Vitis lambrusca cv. Concord (VLRSgt), which pro- 
duced glycosides of stilbenes and glucose esters of 
HCyAs in vitro [54]. Conversely, cinnamyl alcohol de- 
hydrogenase (CAD, Figure 5 '13'), which is involved in 
lignin biosynthesis, displayed quite an opposite profile. No 
other enzymes acting at the branching points of the flavon- 
oid pathway were detected. Together all these results pro- 
vide strong evidence that in ^vQ-veraison stages, carbon 
flows mainly from glycolysis precursors toward not only the 
production of PAC precursors but, and based on weaker evi- 
dence, also hydroxycinnamic acid and stilbene derivatives. 
Alternatively, the production of other pathway end products, 
such as lignins, flavones, flavonols and ACs, seemed blocked 
in the mesocarp by the absence of the corresponding en- 
zymes. These results are highly consistent with the type of 
Phe-derived secondary metabolites accumulation in the de- 
veloping grape berry [3,55-60]. 

Besides ANR, leucoanthocyanidin reductase (LAR) also 
produces PAC precursors; both were expressed during PACs 
accumulation in grapes and were regulated in a temporal- 
and tissue-specific manner [59]. Although LAR was not 



detected in the pericarp of green berries, ANR (Figure 5 
12') could play a relevant role in driving flavonoid biosyn- 
thesis in the epicatechin production in the tissue during 
early berry development (see Figure 5). The general trend 
observed for flavonoid-derived compound synthesis was that 
it was down-regulated in the ripening phase. This stems 
from the observation that CHI and CHS were strongly 
down-regulated at veraison, and remained down-regulated 
until full ripening. Of the remaining flavonoid pathway en- 
zymes, only DFR was detected and displayed almost un- 
changed levels from 15 mm to full ripening. In grape berry 
skins, where ACs heavily accumulated at the onset of ripen- 
ing, no new synthesis of precursors from p-coumaric acid 
apparently occurred in the skin, as observed by the decreas- 
ing profiles of the downsteam enzymes [36]. The presence 
of DFR in the mesocarp during ripening would warrant a 
supply of ACs precursors for synthesis in the sldn via 
ANS and UDP-glucose: flavonoid 3-O-glucosyltransferase 
(UFGT), as suggested by Martinez-Esteso et al [36]. 

Although flavonoid compounds were synthesized in 
the cytosol, they accumulated in vacuoles. It has been 



Martinez-Esteso et al. BMC Plant Biology 2013, 13:167 
http://www.bionnedcentral.conn/1471-2229/13/167 



Page 11 of 20 



thought that the mechanism underlying their transport is 
mediated through gluthatione-S -transferase (GST) and 
ABC transporters [61,62], although GST-independent 
transport also occurs (reviewed by Terrier et aL, [63]). 
High-throughput studies during development in skin- 
colored cultivars have revealed the expression of a GST as- 
sociated with AGs accumulation [14,36,64]. Here, seven 
GSTs belonging to the tau, phi and lambda classes were 
found to be deregulated in flesh during berry develop- 
ment. As Additional file 12 depicts, only one tau class, as- 
sociated with AC accumulation, was also detected in berry 
skin, and none coincided with the four tau class isoforms 
associated with stilbenoid accumulation in elicited cell cul- 
tures [65]. The tau and phi GSTs, in addition to flavonoid 
binding, also perform known auxin and cytokinin binding 
activities [66,67], and a role in plant growth and develop- 
ment has been suggested [62]. 

Three isoenzymes of an isoflavone reductase-like protein 
(IFRL) have been reported as being up-regulated during 
ripening in the mesocarp, which confirms previous results 
obtained in the grape berry mesocarp and skin [31,36]. The 
deregulated IFRLs showed a high sequence homology 
with Eugenol synthase2 from Clarkia brevweri (CbEGS2) 
(79%) and with phenylcoumaran benzylic ether reductase 
(PCBER) from Populus trichocarpa (PtPCBER) (82%). En- 
zymatic activities reported for CbEGS2 show a preference 
on coniferyl acetate over dehydrodiconiferyl alcohol (DDC) 
as a substrate [68] which suggests that the sequences cur- 
rently characterized as PCBER enzymes, as well as other 
phylogenetically related sequences, might prefer coniferyl 
acetate or other related substrates than DDC. In fact, in- 
creased metabolite eugenol during ripening, one of the 
volatiles present in grape berries, has been reported in 
berries of cv. Muscat Hamburg [69]. This increase ran 
in parallel to our quantitative protein profiles for IFRLs. 
Further experiments are required to test these hypotheses. 

3. Photosynthesis, respiration and fermentation 

Apart from the photosynthates which translocate from 
leaves, grape berries are capable of local photosynthetic 
assimilation. A large number of them was detected and 
quantified through these iTRAQ experiments, thus pro- 
viding quite a complete picture of the grape berry photo- 
synthetic machinery during berry development at the 
protein level for the first time (Additional file 13). Almost 
every enzyme involved in the Calvin-Benson cycle was 
identified and quantified. They displayed decreasing abun- 
dance profiles throughout the berry formation and ripen- 
ing processes, with a sharper decrease from 7-to-15 mm 
or 15 mm-to-VlOO. Moreover, a broad representation of 
the light reactions proteins was also detected and quanti- 
fied (Additional file 13). The general profile trend for these 
proteins differed for those involved in Calvin-Benson dur- 
ing berry formation, with light reaction proteins displaying 



a moderate increase in abundance during the first growth 
period which leveling off from 7 mm-to-15 mm. Abun- 
dance reduced during fruit ripening in the same way it did 
for the carbon fixation proteins. These results are consist- 
ent with studies reported at the physiological, transcrip- 
tional and enzyme activity levels, which indicated that the 
capability of grapes to perform photosynthesis diminished 
before veraison and during the ripening process [47,70,71]. 
Likewise, they can correlate largely with transcriptomic 
studies [14,15,17]. The phasing-out of the Calvin-Benson 
cycle in relation to light reaction proteins during the first 
growth period suggests that an important part of the energy 
and reducing power can be used for processes other than 
carbon assimilation, such as the biosynthesis of aromatic 
amino acids and flavonoids, synthesis of proteins or cell 
division, which are highly active during berry formation. 
Indeed, grape berries cannot achieve a net carbon assimi- 
lation rate through photosynthesis at any development 
stage, while photosynthesis would serve to balance the 
loss of CO2 by respiration [72]. 

As noted in the DIGE experiments, oxygen-evolving 
enhancer (OEE) proteins were up-regulated from 7 mm-to- 
15 mm [31]. Subunits of 33 KDa (OEE33) and 16 KDa 
(OEE3) increased more markedly than most PSI, PSII, 
LHCI and LHCII proteins. As hypothesized in a previous 
proteomic study [31], an increase in the OEE compo- 
nents could respond to a redirection of unused electrons 
toward O2 modulating H2O2 formation, which is a 
widely accepted theory when NADP^ is in short supply 
[73,74] and may occur at the end of the berry formation 
phase when many of the aforementioned biosynthetic 
processes slow down, or even stop. Interestingly, two 
fibrillin isoforms (see Additional file 11, 'Other proteins 
of interest') displayed that the protein profiles ran in pa- 
rallel to those of photosynthetic machinery proteins (see 
Additional file 13) with a 1.5- and 2-fold increase for 
each one from 7 mm-to-15 mm, which decreased after 
veraison. Fibrillins are lipid-binding proteins involved in 
the ABA-mediated photoprotection of PSII [75]. The 
fibrillin profile in the grape berry supports the PSII 
stabilization theory for the redirection of electrons and 
H2O2 formation control when it diminishes carbon reduc- 
tion dramatically through the Calvin cycle. In line with 
these findings, antioxidant chloroplastic enzymes were 
simultaneously enhanced (see Additional file 11, 'Stress'). 

As regards respiratory activity, several proteins involved 
in the mitochondrial respiratory chain and oxidative phos- 
phorylation have been detected to be deregulated during 
development (Additional file 9C). The profiles for these 
proteins, along with several subunits for the ATP synthase 
and Complex I NADHiUQ reductase, increased during 
the first growth period, but decreased from the onset of 
ripening, which agrees with the respiratory activities de- 
scribed during development [71,72]. The adenylate and 
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phosphate transporters also displayed similar protein pro- 
files from the onset of ripening (see Additional file 9A). 
Other proteins of the respiratory chain from Complex III 
and cytochrome c remained unchanged or only increased 
slightly. Profiles were consistent with progressive electron 
transport deceleration during ripening. When mitochon- 
drial respiration is compromised, the ATP synthase is able 
to reverse and consume ATP, which serves to maintain 
the mitochondrial membrane potential. This activity can 
deplete ATP and precipitate cell death. This process can 
be mitigated by mitochondrial protein IFi, an endogenous 
ATP synthase inhibitor [76]. An ATPase inhibitor protein 



was found to be strongly up-regulated from 7 mm- 
to-15 mm and during ripening. Here the results also reveal 
that the protein levels of NAD-ME (mME, Figure 6 '24', 
see group 4 'Carbohydrate metabolism') peaked at pre- 
veraison, but declined during ripening. mME was able 
to supply high levels of NADH to the mitochondrial 
matrix when malate was abundant. This reducing environ- 
ment could activate non-phosphorylating dehydrogenases 
(reviewed by Sweetman et al [77]) when the phosphorylat- 
ing pathway is limited, which coincides with the oxidative 
burst observed at veraison [78]. Instead of the alternative 
oxidase as a bypass of proton translocating Complex-I, a 
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Figure 6 Metabolic pathways involved in malate metabolism in grape fruits during development. The protein levels of deregulated 
enzymes are shown in coloured squares, indicating the change of expression (log2 ratio) for each developmental stage in relation to the 15 mm 
stage. In sequence order (left to right), stages are displayed from FS, 4 mm, 7 mm, 15 mm, V-100, 1 10 g/l, and 140 g/l. Different isoforms or 
subunits of proteins are shown as different rows. The enzymes names for each biochemical step are indicated in red if differentially expressed 
and identified, or in black if not detected. Black solid and dashed arrows indicate enzymatic reactions and metabolite movement, respectively. 
SuSy, sucrose synthase; SuPSy, sucrose phosphate synthase; UDP-GP, UDP- glucose pyrophosphorylase; PFK-1: Phosphofructokinase-1; PFP, 
pyrophosphate-dependent phosphofructokinase; FK, fructokinase; HK, hexokinase; PGM, phosphoglucomutase; FBPA, fructose bisphosphate 
aldolase; TPI, triose phosphate isomerase; G3PDH, glyceraldehide-3-phosphate dehydrogenase; PGK, phosphoglycerate kinase; ENO, enolase; PGM, 
phosphoglycerate mutase; PK, pyruvate kinase; PEPC, phosphoenolpyruvate carboxylase; PEPCK, phosphoenolpyruvate carboxykinase; PPDK, 
pyruvate orthophosphate dikinase; cMDH, cytosolic malate dehydrogenase; cME, cytosolic malic enzyme; PDC, pyruvate decarboxylase; ADH, 
alcohol dehydrogenase; vINV, vacuolar acid invertase; VvHT6, hexose transporter 6; PDH, pyruvate dehydrogenase; mMDH, malate dehydrogenase 
mitochondrial; mME, mitochondrial malic enzyme; FH, fumarate hydratase; ATP-CL, ATP-citrate lyase; GS, glutamine synthase; GDH, glutamate 
dehydrogenase; DTC, dicarboxylate-tricarboxylate carrier. 
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mitochondrial uncoupling protein (MUP) was high at pre- 
veraison and down-regulated after the onset of ripening, 
which could contribute to the re-oxidization of coenzymes 
irrespectively of the energy status of berry cells. 

These data provide further evidence for the decreasing 
respiration rates of the mitochondria as ripening proceeds. 
In such a situation, the TCA cycle is expected to become 
progressively inhibited by respiratory control, thus trigger- 
ing the accumulation of malate in the cytosol. Conse- 
quently, part of the malate would be diverted toward 
ethanol fermentation, as supported by pyruvate decarb- 
oxylase accumulation (PDC, Figure 6 '21', see group 4 
'Carbohydrate metabolism') and alcohol dehydrogenase 2 
(ADH2) during ripening, (Figure 6 '22b', see Group 4 
'Carbohydrate metabolism'). ADH2 is the isoform respon- 
sible for ethanol fermentation at ripening [79]. In straw- 
berry, ethanol fermentation can occur aerobically in 
ripening fruit if acidity increases due to, for instance, lim- 
ited respiratory activity, with malate as a putative carbon 
source [80]. The above situation may lead to a shift in the 
cytosolic NAD(P)H:NAD(P)^ balance to a reduced form. 
Interestingly, NADP-dependent FMN quinone reductase 
(BQR) has been reported to markedly increase by up to 4- 
fold from V-lOO-to-110 g/1 and to maintain its levels dur- 
ing ripening. BQR enzymes have not been characterized in 
grape berries and only poorly so in plants [81]. BQRs allow 
for the two-electron reduction of quinones to the hydro- 
quinone form to avoid the generation of one-electron re- 
duced semiquinone, which is known to cause oxidative 
stress [82], while an isoform in yeast has been shown to act 
as a NAD (P) -redox sensor for oxidative stress [83]. Hence, 
BQR becomes an interesting protein to study its putative 
implication in complex redox regulation in the ripening 
grape berry. In addition, BQR has been shown to be one of 
the most abundant proteins in the grape berry mesocarp 
[31]. 

4. Carbohydrate and malate metabolism 

Efforts have been made in recent decades to explain the 
metabolic enzymes controlling the sugar/acid balance 
during grape development, a major quality trait of grape 
berry flesh. Sucrose produced through photosynthesis in 
the mesophyll of mature leaves is loaded into the 
phloem and unloaded into sink berries throughout their 
development and ripening stages (reviewed by Boss and 
Davies, [7]). From fruit set to veraison, malate is stored 
in vacuoles as a major end product of imported sucrose; 
thereafter malate is consumed, while large amounts of 
sucrose-derived glucose and fructose are stored in the vacu- 
ole. As demonstrated by Zhang et al [50], the unloading 
pathway shifts from symplastic to apoplastic at the onset of 
ripening, and is accompanied by a concominant increase in 
the expression and activity of cell wall invertases (cwINV). 
Proteomic data [31] have provided evidence for a switch 



between sucrose synthase (SuSy) and vacuolar acid invertase 
(vINV) before veraison, thus supporting a supposed ftinc- 
tional mechanism between SuSy and INV for unloading 
sugars during grape berry development and ripening, re- 
spectively [84], as demonstrated in tomato [85-87]. In the 
current experiments, cwINV was not identified, but a vINV 
(Figure 6, '23') was strongly up-regulated (14-fold) from 7- 
to-15 mm and protein levels were maintained during ripen- 
ing. The transcripts and protein levels of two vINV 
(VvGINl and VvGIN2) accumulated before veraison and 
decreased during fruit ripening [4,17], although vINV activ- 
ity peaked before veraison and remained constant on a per 
berry basis [4]. Moreover, the natural depletion of vINV in 
the Steuven grapevine hybrid contributed to increase su- 
crose in the maturing berry [88]. The above cumulative evi- 
dence indicates that, apart from the important role of 
cwINV during ripening [50], vINV also plays a relevant role 
in driving the import of sugars during ripening. Further- 
more, SuSy (Figure 6, '1') has been identified as two 
isoforms. One was abundant since the beginning of develop- 
ment to 7 mm, but then decreased abruptly before veraison, 
displaying the opposite trend to that of vINV (Figure 6, 'la'). 
The second SuSy isoform (Figure 6 'lb') appeared to be up- 
regulated at the onset of ripening, while the levels were 
maintained thereafter until ftill ripening. Since the SuSy ex- 
pression was induced by a sugar-sensing mechanism [89], 
its decrease before veraison suggests a lower sucrose con- 
centration in the cytoplasm of sink cells, probably due to 
plasmodesmata blocking taking place at the shift of the 
unloading pathway [50]. The first isoform, which was also 
detected in a previous study [31], is a clear candidate to be 
involved in both sugar unloading and metabolism during 
the first growth period, while the second isoform would be 
responsible for the function during ripening, together with 
vINV. The present iTRAQ experiments revealed a third en- 
zyme involved in sucrose unloading and sink strength, su- 
crose phosphate synthase (SuPSy) (Figure 6 '2'), which was 
up-regulated almost 4-fold at the onset of veraison and 
remained at constant levels during ripening. The joint action 
of SuSy, SuPSy and INV would result in dynamic sucrose 
synthesis and degradation cycles, referred to as 'ftitile cycles,' 
which occurred in the plant cells regulating flow intensity 
and direction. Sucrose futile cycles are key mechanisms for 
unloading and storing sugars into ripening tomato fruit [87] 
that, as in grapes, stops the symplastic unloading pathway at 
the onset of ripening [50]. 

Malate is the main organic acid in grape berries and 
shows an accumulation pattern, which peaks in the pre- 
veraison stage before ripening is triggered [90]. Its metab- 
olism operates tightly in parallel with sugar fate and 
carbohydrate metabolism in grape berry cells in the meso- 
carp. The abundance profiles for the enzymes involved in 
carbohydrate metabolism tie in with those found in a pre- 
vious DIGE proteomic study [31]. Moreover, this iTRAQ 
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large-scale study has allowed to increase the coverage of 
the metabolic pathways driving the sugar/acid balance 
during development, and more importantly, also of those 
enzymes controlling the glycolysis/gluconeogenesis and 
synthesis/degradation of malate. A previous study has 
highlighted the relevance of pyrophosphate-dependent 
phosphofructokinase (PFP, Figure 6 '8') as the catalyst of 
the phosphorylation of fructose-6 phosphate during green 
development. As discussed before (Figure 4), the low levels 
of a sPPase and the occurrence of PPi releasing reactions 
during green development support the availability of PPi 
for PFP. PFP, unlike phosphofructokinase- 1 (PFK-1, Figure 6 
7'), is not negatively regulated by end products such as PEP 
or ATP, thus allowing for the glycolytic flow to proceed 
driven simply by substrate availability. Here the profile of 
PFP reveals that, after dropping just before veraison, the 
levels are no longer recovered. PFK-1, detected in the ripen- 
ing phase only, was strongly up-regulated at veraison and 
its levels were held until the end of ripening. Hence, during 
ripening glycolysis would be controlled strictly by accumu- 
lated PEP and by the energy status of mesocarp cells. 

The next key event was the usage of PEP during green 
development. Several lines of evidence support two major 
uses: conversion into malate via oxalacetate (OAA) for 
storage in vacuoles [17,69,91] and respiration [72]. The 
profile of the corresponding enzymes acting on PEP indi- 
cated that both pathways are open (Figure 6). 

The two enzymes converting PEP into OAA, PEP carb- 
oxylase (PEPC, Figure 6 '20') irreversibly and PEP 
carboxykinase (PEPCK, Figure 6 15') reversibly, were 
detected here, but with quite different profiles. One isoform 
of PEPC was abundant during early development, but 
dropped before veraison, while another isoform displayed 
an abundance peak at veraison. Such profiles are consistent 
with those for gene expression [14,77,92] and enzyme activ- 
ity [93,94] in developing grape berries, thus providing evi- 
dence for its role in malate synthesis at the protein level. 
Instead, the abundance of PEPCK isoforms started to in- 
crease in the last green development stages, and their level 
was held until the end of ripening. Despite the detection of 
transcripts and enzyme activity in ^VQ-veraison berries 
[14,77,95], the profiles found herein do not support a key 
role of PEPCK in malate synthesis, as previously suggested 
[96], rather they would be important during ripening, as 
discussed below. The subsequent conversion of OAA into 
malate in the green stages by cytosolic malate dehydrogen- 
ase (cMDH) is also well supported by the cMDH profile 
(Figure 6 19'), by the active glycolysis and PEPC in this 
phase as continuous sources of NADH and OAA, respect- 
ively, and by the thermodynamics of the reaction itself. Al- 
though no individual measurements of cMDH activity have 
been reported, the profiles found herein, as compared to 
those of mitochondrial MDH (mMDH), suggest that it may 
contribute largely to the high levels of total MDH activity 



reported during early development [97]. Taken together, 
the present results support the notion that a major carbon 
flow in malate synthesis in the green phase of grape berry 
development occurs via PEPC and cMDH until its max- 
imum accumulation at ^vQ-veraison, 

In respiratory usage, PEP is directly converted into pyru- 
vate through pyruvate kinase (PK, Figure 6 16') and is 
transported to mitochondria for complete oxidation. The 
respiratory use of PEP is supported by the PK profile, which 
is more abundant during early development than at pre- 
veraison, when it is minimal. The profiles of pyruvate de- 
hydrogenase (PDH, Figure 6 '27') and MDHm in that phase 
are flat. The CO2 evolution profiles in grape berries, higher 
early on in development to decrease at the beginning of rip- 
ening [72], correlate better with glycolytic enzymes, includ- 
ing PK, than with mitochondrial enzymes PDH and 
mMDH, which suggests that the respiratory machinery is 
not the rate-limiting step. Since both glycolysis and photo- 
synthesis are operative, the carbon flow through the former 
and the intensity of the latter (see the profiles above) may 
be major determinant factors of CO2 evolution. 

Pyruvate might potentially be converted into malate by 
the cytosolic NADP-dependent malic enzyme (cME, 
Figure 6 18'). Yet despite cME catalyzing a reversible reac- 
tion, malate decarboxylation is thermodynamically favored. 
So this enzyme is considered to be involved in malate deg- 
radation during ripening [98]. cME levels rose halfway 
through green development, while malate accumulated and 
dropped at mid-ripening. In this sense, cME may not be 
assigned a unique role in malate degradation during ripen- 
ing, when large amounts of this compound are released 
from the vacuole. One such role can be speculated to form 
part of a biochemical pH-stat, which also involves MDH, 
PEPC and PEPCK [99]. 

At the inception of ripening, net malate accumulation 
switches to net degradation [8] due to its release from 
the vacuole and its use in different pathways. Gluconeo- 
genesis has been suggested to occur in grapes in ripen- 
ing stages [100]. As mentioned above, cMDH abundance 
decreased notably at ^vQ-veraison, while two PEPCK 
isoforms showed a moderate accumulation from the 
^vQ'Veraison stage. The increase in the PEPCK tran- 
scripts [14,77], enzyme activity during ripening [95] and 
the proteins found herein supports its role in gluconeo- 
genesis. However another pathway involving cME and 
pyruvate diphosphate dikinase (PPDK, Figure 6 17') may 
also operate in gluconeogenesis in plants. A transcript 
encoding a putative PPDK has been seen to increase 
throughout berry development, while no changes in the 
expression of the ME isoforms has been detected [77]. 
The accumulation profiles of cME found herein, from 
mid-green development to mid-ripening, and PPDK, 
which peaked at veraison, are consistent with this role in 
a transient gluconeogenic phase at the beginning of 
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ripening. If we bear in mind the switch of the sugar 
unloading pathway at pre-veraison [50], it can be hy- 
pothesized that gluconeogenesis at early ripening stages 
is required to compensate the transient decay of sugars 
from the time that plasmodesmata are blocked to the 
full operation of the transporter-mediated apoplastic 
sugar unloading. The accumulation of PFK-1 and down- 
stream glycolytic enzymes, including PK, from veraison 
to full ripening indicates that the pathway may flow to- 
ward pyruvate once the sugar import is restored once 
again. Early biochemical and physiological studies pro- 
vided evidence that gluconeogenesis occurs, but that it is 
not a major pathway of malate degradation throughout 
whole berry ripening [90]. Accordingly, the profiles of 
cMDH, PEPCK, cME and PPDK support gluconeogene- 
sis in early, but not late, ripening stages. 

Respiration during ripening is strongly supported by the 
corresponding profiles of the TCA cycle enzymes and, at 
least during the first half of the ripening phase, by the pro- 
files of respiratory complexes (see group 3 above about 
'Respiration'). Large quantities of malate released from 
vacuoles from veraison and pyruvate, obtained from either 
malate or glycolysis, may be transported directly to mito- 
chondria to feed the TCA cycle, to produce ATP and to 
maintain the respiratory flux in fruit cells. The up- 
regulation of mMDH (Figure 6 '25') ties in with increasing 
protein levels and activity at post-v eraison for MDH 
[97,101]. The use of pyruvate as a respiratory substrate is 
supported by the increasing PDH profile (Figure 6 '27'). 
This scenario is in agreement with our hypothesis for a 
change in the malate degradation metabolism from being 
supported mainly by respiration at early ripening stages to 
then occur by fermentation in the cytosol at late ripening. 
In the late ripening stages, the respiration rate may de- 
crease and induce ethanol synthesis (see group 3 above 
about 'Fermentation'). 

Conclusion 

Applying iTRAQ to study grape berry development has 
allowed the identification and quantitation of 411 and 
630 proteins in the green and ripe growing phases, re- 
spectively. These longer lists of proteins complement a 
previous gel-based proteomic study with a better prote- 
ome coverage and they better connect the two growing 
phases by analyzing a common time point in develop- 
ment, particularly 15 mm. This technique allowed the 
detection of another key point in development, 15 mm- 
to-VlOO, where most of the dramatic changes at the 
protein level occurred. The obtained results led to a com- 
prehensive study of grape berry development that sup- 
ports and complements a previous proteomic analysis 
[31]. This large-scale proteomic study as a hypothesis-free 
approach provides quite a complete view of the major 
and important pathways which evolve during fruit 



development, thus providing a better understanding of 
berry development and ripening physiology. These find- 
ings help provide an understanding of the metabolism and 
storage of sugars and malate, energy-related pathways 
such as respiration, photosynthesis and fermentation, and 
the synthesis of polyphenolics as major secondary metabo- 
lites in grape berries. They all largely determine final berry 
quality and are of paramount importance for the viticul- 
tural industry. A similar approach at the protein level is 
now feasible to study the effect of other variables of inter- 
est (environmental factors and cultural practices) on these 
pathways. Alternatively, the key steps identified in this 
study, such as the PFP-PFK or SuSy-INV switches, among 
others, can be targeted under multiple conditions to finely 
characterize their influence on the final sugar/acid balance 
in ripe fruit. Finally, some proteins underwent major 
changes at specific developmental stages; thus, they can be 
used as novel protein biomarkers of berry development, 
which have not been detected to date. Consequently, this 
may help open up new lines to explore the parameters 
controlling grape berry development and ripening. 

Methods 

Plant material 

Grape berries {Vitis vinifera L. cv. Muscat Hamburg) 
were collected from the experimental vineyard at the 
Instituto Murciano de Investigacion y Desarrollo 
Agrario (Torrepacheco, Murcia, Spain) in 2006. Berries 
were sampled from four selected vines and were col- 
lected at seven different developmental stages from 
fruit set until full ripening during both growing sea- 
sons. Individual grapes were developmentally staged 
according to the different berry growth pattern param- 
eters. Green berries were classified into fruit set (stage 
FS) and then according to the equatorial diameter of fruit 
in 4, 7 and 15 mm. The beginning of the second berry 
growth phase, labeled 100% veraison (stage V-100), was 
assessed visually as 100% berry surfaces turned pink in. 
Ripening berries were classified according to their esti- 
mated density by flotation in different NaCl solutions: 
stage 110 was berries that sink in 110 g/1, but float in 120 
g/1, while stage 140 was berries that sink in 140 g/1. Since 
berries in a bunch do not develop uniformly, berries were 
harvested from different bunches in the same plant on 
each sampling day and were then sorted according to their 
developmental stage, considered to be a sample. Four se- 
lected vines were sampled during berry development and 
ripening to make four biological replicates per stage. Sam- 
pled berries were immediately frozen in liquid nitrogen 
after detaching, except for the 110- and 140-staged berries, 
which were classified previously according to density. All 
samples were stored at -80°C until use. A parallel set of 
sampled berries was refrigerated after detaching and was 
transported to the laboratory to determine color index. 
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juice pH, total acidity and °Brix, which showed typical pro- 
files for grapevine berries. The methods and results are 
reported elsewhere [69]. 

Protein extraction 

The three main berry tissues, seed or endocarp, flesh or 
mesocarp and skin or exocarp, were differentiated at the 
different developmental stages, and their dissection was 
only possible after certain stages. Prior to total protein 
extraction, seeds were removed from berries from stage 
4 mm onward, while exocarp tissue was peeled away in 
berries from stages V-100 onward. Thus, green berry 
proteins and ripening berries proteins were extracted 
from the pericarp and the mesocarp, respectively. The 
tissue of a pool of berries from the same vine was 
ground to a fine powder in a mortar with liquid nitrogen 
and 4 g of ground tissue were used to prepare a protein 
extract. Protein extracts were obtained as described in 
Martinez-Esteso et aL [31]. 

The protein was quantified by the Bradford dye- 
binding method [102] with bovine serum albumin used 
as a standard and an equal amount of protein from each 
staged replicate was pooled for isobaric labeling. 

Isobaric peptide labeling 

For each developmental stage, a volume corresponding to 
100 [ig of protein was precipitated with 10 volumes of acet- 
one at -20°C overnight. After centrifugation for 10 min at 
15300 X the protein pellet was dissolved in 60 [A of 
iTRAQ dissolution buffer (Applied Biosystems) containing 
0.2% (w/v) SDS. Proteins were reduced in 3 mM tris- 
(2-carboxyethyl) phosphine (TCEP) and were incubated for 
1 h at 60°C. After cooling samples to RT, cysteine residues 
were blocked with 2 [A 200 mM methylmethanethiosulfate 
(MMTS) by incubating at RT for 10 min. Samples were di- 
luted by adding 250 [A of iTRAQ dissolution buffer. Then, 
10 [ig of proteomics grade modified trypsin (Sigma), 
dissolved in the same buffer, were added to each vial and di- 
gestion was allowed to proceed overnight at 37°C. After- 
ward, a small pellet remained and 5 [ig of proteomics grade 
modified trypsin (Sigma) dissolved in iTRAQ dissolution 
buffer were added to each vial and allowed to digest at 
37°C for 3 h. The resulting tryptic peptides were 
vacuum-concentrated and re-suspended in 30 [A of iTRAQ 
dissolution buffer. The labeling reactions were done fol- 
lowing the manufacturers recommendations by adding 
one iTRAQ reagent, 114.11123, 115.10826, 116.11162 or 
117.11497, previously dissolved in 70 [A of pure ethanol, to 
each protein sample vial. The labeling reaction was stopped 
after 1 h of incubation at RT by adding 1 ml of a buffer 
containing 10 mM K2HPO4 and 25% acetonitrile (ACN), 
pH 2.7, to each vial. In one iTRAQ experiment, reagents 
were used to label the FS stages 4 mm, 7 mm and 15 mm, 
and to label the stages 15 mm, V-100, 110 g/1 and 140 g/1 



in another. Then, the four samples were pooled and ad- 
justed to pH 3.0 with concentrated phosphoric acid. As all 
these experiments were carried out in a different laboratory, 
the downstream workflow procedures differed slightly. 

Peptide fractionation by strong cation exchange 

The pool of labeled samples was fractionated by strong 
cation exchange chromatography (SCX). Samples were 
separated using an Akta Purifier (GE Healthcare) medium 
pressure liquid chromatography system equipped with a 
Mono S PC, 1.6 mm x 50 mm column (GE Healthcare) at 
a flow rate of 0.1 ml/min (green stages), or a BioCAD 
workstation (Applied Biosystems) utilizing a 100 x 4.6 
mm polysulfoethyl aspartamide column (PolyLC Inc, 
Columbia, MD, USA) at a flow rate of 0.5 ml/min (ripe 
stages). First, the sample was loaded into the column, 5- 
fold diluted in buffer A (10 mM K2HP04/25% ACN pH 
2.7). Afterward, the column was washed with buffer A for 
20 min and peptides were eluted with a two-step gradient: 
first a linear gradient of 5-35% buffer B (0.5 M KCl in 10 
mM K2HP04/25% ACN pH 2.7) for 30 min, followed by a 
linear gradient of 35-100% buffer B for 60 min. The elu- 
tion of peptides was monitored at 280 nm and fractioned 
into 0.1 ml (green stages) or 214 nm and 0.5 ml (ripe 
stages) throughout the chromatographic run. For the 
green stages, the 77 collected fractions were further re- 
duced to 21 by pooling sets of three or four consecutive 
fractions. The 21 resulting fractions were vacuum- 
concentrated, resuspended in 100 [A 5% ACN/0.5% 
trifluoroacetic acid (TEA) and desalted with PepClean™ 
C-18 Spin Columns (Thermo Fisher Scientific, Rockford, 
IL) according to the manufacturers recommendations. 
The peptides eluted in 40 \A 70% ACN/0.1% formic acid 
(FA) were dried under vacuum and re-suspended in 18 [A 

0. 1% FA. For the ripe stages, 61 fractions were collected, 
but only 18 containing the eluted labeled peptides mea- 
sured by optical density monitoring at 214 nm were chosen 
for the analysis in a 2-hour LC-MS/MS program. The frac- 
tionated samples were reduced to 150 [A in a speed-vac 
(Thermo-Savant, Holbrook, NY, USA) and were transferred 
to autosampler tubes (LC Packings, Amsterdam, The 
Netherlands). 

Reverse phase chromatography 

An integrated system consisting of a Famos Autosampler, a 
Switchos switching pump and a UltiMate micropump (LC 
Packings, Amsterdam, Netherlands) was used for the re- 
verse phase chromatography on the selected SCX fractions. 
For the green stages, 5 [i\ of the tryptic peptides were 
pre-concentrated in a CI 8 PepMap guard column (300 [xm 

1. d. X 5 mm 5 (im, 100 A, LC Packings, Amsterdam, The 
Netherlands) at 40 (il/min for 3 min in 0.1% FA, followed 
by elution in a C18 PepMap (75 [im i.d. x 15 cm, 3 (im, 
100 A, LC Packings, Amsterdam, The Netherlands) using a 
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120 min linear gradient from 15 to 50% solvent B; solvent 
A was 0.1% FA in water and solvent B was 0.1% FA in 95% 
ACN. For the ripe stages, one fifth of each SCX fraction 
was desalted in a CI 8 PepMap guard column (300 [im 
i.d. X 5 mm 5 (im, 100 A, LC Packings, Amsterdam, 
The Netherlands) at 50 (il/min for 15 minutes. HPLC 
buffers consisted of Buffer A -2% ACN, 0.1% FA and 
Buffer B -98% ACN and 0.1% FA. Peptides were sepa- 
rated ina manually packed 75 [im x 15 cm CI 8 column 
(Magic C18Aq, 5 (im, 100 A, Michrom Bioresources 
Inc., Auburn CA, USA) using an 85 min gradient of 
5%-75% buffer B flowing at 250 nl/min for fractions 
37-48 and a 35 min gradient of 5%-75% buffer B. 

Mass spectrometry 

The eluent was sprayed directly into either a QSTAR XL 
System (green stages) or a QSTAR Pulsar I (ripe stages) 
mass spectrometer equipped with a nanospray source 
(Applied Biosystems/MDS SCIEX Concord, ON Canada). 
The QSTAR operating software Analyst QS vl.l employed 
an information dependent acquisition (IDA) method for 
optimized MS/MS spectra acquisition over a 6-second 
cycle, which was repeated throughout gradient duration. 
The MS-TOF survey scan lasted 1 second over the range 
of 400-1200 m/z targeting ions of charge state 2-4+, 
which exceeded a threshold of 20 counts. The former tar- 
get ions within 100 ppm were excluded for the next 180 
seconds. Each product ion scan lasted 2.5 seconds over a 
range of 100-1500 m/z. Enhance all was turned on for the 
product ion scans. 

Database search and protein quantitation 

Raw data files were processed using Protein Pilot vl.O 
with the Paragon Search and ProGroup Algorithms™ 
(Applied Biosystems/MDS Sciex Foster City, CA USA) 
for both tryptic peptide identification and quantitation. 
The peptides and corresponding relative abundances 
were obtained in ProteinPilot using a confidence cutoff 
(called a 'Prot Score') of >1.0 (>90%) and >1.3 (>95%) for 
the experiments of the green and ripe stages, respect- 
ively. Database searching for each sample was done 
against the NCBInr protein database without taxonom- 
ical restrictions, trypsin, MMTS as a fixed modification 
and the iTRAQ label as a variable modification. Only 
the proteins identified with at least 2 different peptides 
and p<0.05, and quantified with a ratio of >1.5 and 
p<0.05, were considered. The former p-value related to 
the protein score cutoff in the identification, while the 
latter p-value related to the iTRAQ ratio for each quan- 
tified protein and was computed from the ProGroup 
Algorithm in the ProteinPilot software as a measure of 
its statistical significance. The sequences the from kera- 
tins, trypsin and species other than plants were not 
considered. 



Bioinformatic functional analysis 

Since the translated ORFs from the grapevine genome pro- 
jects [34,35] have been deployed in public databases with- 
out descriptions and annotations to date, searches often 
match an undescribed amino acid sequence. The Blast2GO 
v2.4.0 application [43] has been used to automatically as- 
sign protein description and take-up annotations from 
homologous sequences of public databases, which have 
then been manually reviewed and enriched, if possible. A 
file of the FASTA format sequences of the identified and/or 
quantified protein set was batch-retrieved from the NCBI 
website. Blast2GO was fed with the FASTA file and was 
run to: first incorporate the sequence description by 
performing a BLASTp search against NCBInr (e-value cut- 
off of 1 x 10'^^, 20 for the retrieved number of BLAST hits, 
33 for the HSP (highest scoring pair) length cutoff); second 
to map the GO, EC and Interpro terms; then to annotate 
the sequences (E- Value Hit-Filter of 1 x 10'^, a Hsp-Hit 
Coverage Cutoff of 0, an Annotation Cutoff of 55, and a 
GO Weight of 5). The automatic annotation performed by 
Blast2GO was manually revised to guarantee accurate as- 
signment. Fisher exact tests [45] were performed to find 
the significantly enriched GO terms (FDR <5%). Briefly, the 
frequencies of the annotation terms were compared be- 
tween a reference list containing all the quantified se- 
quences in the experiment and a subset list of sequences 
selected through an arbitrary quantitation ratio cutoff of 
1.5-fold (ratio <0.6 for the down-regulated and >1.5 for the 
up-regulated proteins). Subsets were selected at each devel- 
opmental stage transition analyzed: FS-to-4 mm, 4 mm- 
to-7 mm, 7 mm-to 15 mm, 15 mm-to-VlOO, VlOO-to-110 
and llO-to-140. 
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Additional file 1: List of proteins identified and quantified for the 
iTRAQ experiment of the green stages, including the protein and 
peptide output files. 

Additional file 2: Processed output file (Additional file 1), including the 
protein description annotation by Blast2G0 v2.4.0 application [43]. 

Additional file 3: List of proteins identified and quantified for the 
iTRAQ experiment of the ripe stages, including the protein and 
peptide output files. 

Additional file 4: Processed output file (Additional file 3), including the 
protein description annotation by Blast2G0 v2.4.0 application [43]. 

Additional file 5: Distribution of the annotations of biological 
processes and cellular components. 

Additional file 6: File containing the GO-terms annotated by the 
Blast2G0 V2.4.0 application [42] for the proteins identified in the 
iTRAQ experiment of the green stages. 

Additional file 7: File containing the GO-terms annotated by 
Blast2G0 V2.4.0 application [42] for the proteins identified in the 
iTRAQ experiment of the ripe stages. 

Additional file 8: The GO terms enriched and sequences annotated 
with such terms in the up- or down-regulated subsets of proteins in 
each grape berry developmental stage step analyzed. 
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Additional file 9: Profiles of protein functional clusters during berry 
development [103]. 

Additional file 10: Methione synthesis and connection with folate 
metabolism in grape berries during development. 

Additional file 11: Results and discussion about the functional 
groups 'Nitrogen and amino acid metabolism', Terpenoid 
metabolism', 'Signaling and hormone', 'Stress', 'Protein synthesis', 
'Protein degradation', 'Protein processing', 'Cell division and 
growth, biogenesis', 'Defense proteins', 'Other proteins of interest'. 

Additional file 12: Phylogram of the grapevine GSTs identified in 
grape berry during development [104]. 

Additional file 13: Scenario of the photosynthetic machinery of 
both the light and dark phases during grape berry development. 



Abbreviations 

ACs: Anthocyanins; DIGE: Difference gel electrophoresis; FS: Fruit set; 
VI 00: Full color char^ge at veraison; GO: Gene ontology; GST: Gluthatione- 
S-transferase; HSPs: Heat shock proteins; iTRAQ: Isobaric tags for relative and 
absolute protein quantitation; HPLC-MS/MS: Liquid chromatography coupled 
with tandem mass spectrometry; cMDH: Malate dehydrogenase cytosolic; 
mMDH: Malate dehydrogenase mitochondrial; PFK-1: Phosphofructokinase 1; 
PACs: Proanthocyanidins; PFP: Pyrophosphate-dependent 
phosphofructokinase; sHSP: Small HSP; SCX: Strong cation exchange 
chromatography; SuSy: Sucrose synthase; TFA: Trifluoroacetic acid; 
2-DE: Two-dimensional gel electrophoresis; UDP-GP: UDP-glucose 
pyrophosphorylase; vINV: Vacuolar acid invertase. 

Authors' contributions 

MJME contributed to the experimental design, grape sampling, protein extraction, 
labeling and SCX fractionation, bioinformatic analysis and biological and data 
interpretation. contributed to grape sampling, protein extraction and labeling. 
MLV carried out the mass spectrometric analysis for the protein identification and 
quantification in the iTRAQ experiment of the green stages. MAP contributed to 
the experimental design and data interpretation. RBM defined the work objectives 
and technical approach, and contributed to the experimental design, bioinformatic 
analysis and data interpretation. This article is part of Maria Jose Martinez-Esteso's 
PhD Thesis. All authors read and approved the final manuscript. 

Acknowledgements 

MJME acknowledges a grant from the University of Alicante. The authors 
gratefully acknowledge the collaboration of the Instituto Murciano de 
Investigacion y Desarrollo Agrario Alimentario (IMIDA) in providing the grape 
berry material. MJME acknowledges a stay in the UVic Genome BC Proteomics 
Centre, with PhD Christoph Borchers as a facility director, in order to perform 
the iTRAQ experiment of the ripe stages analysis. MJME wish to thank Monica 
Elliott and Derek Smith for experimental assistance during the stay. The mass 
spectrometry analysis of the iTRAQ experiment of the green stages was carried 
out by the Proteomics Unit at the Centro de Investigaciones Principe Felipe, 
directed by Dr. Manuel Sanchez del Pino. The University of Alicante Proteomics 
Facility is a member of the Proteored network. This study was supported by 
Genoma Espaha, the GrapeGen project, as a part of a Genoma Espaha-Genome 
Canada joint collaborative research agreement. 

Author details 

^Grupo de Proteomica y Genomica Funcional de Plantas, Dept. Agroquimica 
y Bioquimica, Facultad de Ciencias, Universidad de Alicante, Apartado 99, 
E-03080 Alicante, Spain. ^Grupo de Peroxidasas Vegetales, Department 
Fisiologia Vegetal, Facultad de Biologia, Universidad de Murcia, Campus de 
Espinardo, E-30100 Murcia, Spain. ^Laboratorio de Proteomica, Centro de 
Investigacion Principe Felipe, Av. Autopista del Saler, 16, 46012 Valencia, 
Spain. 

Received: 25 February 2013 Accepted: 8 October 2013 
Published: 24 October 2013 

References 

1 . FAO:food and agriculture organization of the United Nations; 2007. http:// 
faostat.fao.org/site/567/DesktopDefault.aspx?PagelDy4567#ancor 
(16 December 2010, date last accessed). 



2. Coombe BG: Research on development and ripening of the grape berry. 

Am J Enol Vitic ]992, 43:101-110. 

3. Conde C, da Silva P, Pontes N, Dias ACP, Tavares RM, Sousa MJ, Agasse A, 
Delrot S, Geros H: Biochemical changes throughout grape berry 
development and fruit and wine quality. Food 2007, 1:1-22. 

4. Davies C, Robinson SP: Sugar accumulation in grape berries: cloning of 
two putative vacuolar invertase cDNAs and their expression in 
grapevine tissues. Plant Physiol 1996, 1 1 1(l):275-283. 

5. Kanellis AK, Roubelakis-Angelakis KA: Grape. In Biochemistry of fruit ripening. 
Edited by Seymour G, Taylor J, Tucker G. London: Chapman and Hall; 
1993:189-234. 

6. Ruffner HP, Hawker JS: Control of glycolysis in ripening berries of Vitis 
vinifera. Phytochemistry 1 977, 1 6:1 1 71 -1 1 75. 

7. Boss PK, Davies C: Molecular biology of sugar and anthocyanin 
accumulation in grape berries. In Molecular biology and biotechnology of 
the grapevine. Edited by Roubelakis-Angelakis KA. Dordrecht: The 
Netherlands: Kluwer Academic Publishers; 2001:1-33. 

8. Ribereau-Gayon P, Boidron JN, Terrier A: Aroma of Muscat grape varieties. 
JAgrFood Chem 1975, 23:1042-1047. 

9. Marino NC, Tamames EL, Jares CMG: Contribution to the study of the 
aromatic potential of three Muscat Vitis vinifera varieties: identification 
of new compounds. Food Sci Technol Int 1995, 1:105-1 16. 

1 0. Boss PK, Davies C, Robinson SP: Analysis of the expression of anthocyanin 
pathway genes in developing Vitis vinifera L. cv Shiraz grape berries and 
the implications for pathway regulation. Plant Physiol 1 996, 1 1 1 :1 059-1 066. 

11. Tattersall DB, van Heeswijck R, Hoj PB: Identification and characterization 
of a fruit-specific thaumatin-like protein which accumulates at very low 
levels in conjunction with the onset of sugar accumulation and berry 
softening in V. vinifera. Plant Physiol 1997, 114:759-769. 

1 2. Davies C, Wolf T, Robinson SP: Three putative sucrose transporters are 
differentially expressed in grapevine tissues. Plant Science 1999, 147:93-100. 

13. Nunan KJ, Davies C, Robinson SP, Fincher GB: Expression patterns of cell 
wall-modifying enzymes during grape berry development. Planta 2001, 
214:257-264. 

14. Terrier N, Glissant GJ, Barrieu F, Abbal P, Coutre C, Ageoges A, Atanassova R, 
Leon C, Renaudin JP, Dedaldechamp F, Romieu C, Delrot S, Hamdi S: 
Isogene specific oligo arrays reveal multifaceted changes in gene 
expression during grape berry (Vitis vinifera L.) development. 

Planta 2005, 222:832-847. 

15. Waters DLE, Holton TA, Ablett EM, Lee LS, Henry RJ: cDNA microarray 
analysis of the developing grape {Vitis vinifera cv. Shiraz) berry skin. 
Funct Integr Genomics 2005, 5:40-58. 

16. Waters DLE, Holton TA, Ablett EM, Slade LL, Henry RJ: The ripening wine 
grape berry skin transcriptome. Plant Sci 2006, 171:132-138. 

1 7. Deluc LG, Grimplet J, Wheatley MD, Tllet RL, Quilici D, Osborne C, Schlauch 
KA, Schooley DA, Cushman JC, Cramer GR: Transcriptomic and metabolite 
analyses of Cabernet Sauvignon grape berry development. 

BMC Genomics 2007, 8:429. 

18. Fernandez L, Torregrosa L, Terrier N, Sreekantan L, Grimplet J, Davies C, 
Thomas MR, Romieu C, Ageorges A: Identification of genes associated 
with flesh morphogenesis during grapevine fruit development. Plant Mol 
Biol 2007, 63:307-323. 

19. Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, 
Pezzotti M, Delledonne M: Characterization of transcriptional complexity 
during berry development in Vitis vinifera using RNA-Seq. Plant Physiol 
2010, 152:1787-1795. 

20. Castro AJ, Carapito C, Zorn N, Magne C, Leize E, Van Dorsselaer A, Clement 
C: Proteomic analysis of grapevine (Vitis vinifera L.) tissues subjected to 
herbicide stress. J Exp Bot 2005, 56(421 ):2783-2795. 

21 . Vincent D, Ergul A, Bohlman MC, Tattersall EAR, Tllett RL, Wheatley MD, 
Woolsey R, Quilici DR, Joets J, Schlauch K, Schooley DA, Cushman JC, 
Cramer GR: Proteomic analysis reveals differences between Vitis vinifera 
L. cv. Chardonnay and cv. Cabernet Sauvignon and their responses to 
water deficit and salinity. J Exp Bot 2007, 58:1873-1892. 

22. Jellouli N, Ben Jouira H, Skouri H, Ghorbel A, Gourgouri A, Mliki: A 
Proteomic analysis of Tunisian grapevine cultivar Razegui under salt 
stress. J Plant Physiol 2008, 1 65:471-481 . 

23. Martinez-Esteso MJ, Selles-Marchart S, Vera-Urbina JC, Pedreno MA, Bru-Mart 
Inez R: Changes of defense proteins in the extracellular proteome of 
grapevine {Vitis vinifera cv. Gamay) cell cultures in response to elicitors. 

J Proteomics 2009, 73:331-341. 



Martinez-Esteso et al. BMC Plant Biology 2013, 13:167 
http://www.bionnedcentral.conn/1471-2229/13/167 



Page 1 9 of 20 



24. Sarry JE, Sommerer N, Sauvage FX, Bergoin A, Rossignol M, Albagnac G, 
Romieu C: Grape berry biochemistry revisited upon proteomic analysis of 
the mesocarp. Proteomics 2004, 4:201-215. 

25. Vincent D, Wheatley MD, Cramer GR: Optimization of protein extraction of 
mature grape berry clusters. Electrophoresis 2006, 27:1853-1865. 

26. Deytieux C, Geny L, Lapaillerie D, Claverol S, Bonneu M, Doneche B: 
Proteome analysis of grape skins during ripening. J Exp Bot 2007, 
58:1851-1862. 

27. Giribaldi M, Perugini I, Sauvage FX, Shubert A: Analysis of protein changes 
during grape berry ripening by 2-DE and MALDI-TOF. Proteomics 2007, 
7:3154-3170. 

28. Zhang J, Ma H, Feng J, Zeng L, Wang Z, Chen S: Grape berry plasma 
membrane proteome analysis and its differential expression during 
ripening. J Exp Bot 2008, 59:2979-2990. 

29. Karp NA, Feret R, Rubtsov DH, Lilley KS: Comparison of DIGE and post- 
stained gel electrophoresis with both traditional and SameSpots analysis 
for quantitative proteomics. Proteomics 2008, 8:948-960. 

30. Di Carii M, Zamboni A, Pe ME, Pezzotti M, Lilley KS, Benvenuto E, Desiderio 
A: Two dimensional differential in gel electrophoresis (2D-DIGE) analysis 
of grape berry proteome during postharvest withering. J Proteome Res 
2011, 10(2):429-46. 

31. Martinez-Esteso MJ, Selles-Marchart S, Lijavetzky D, Pedreho MA, Bru-Mart 
Inez R: ADIGE-based quantitative proteomic analysis of grape berry flesh 
development and ripening reveals key events in sugar and organic acid 
metabolism. J Exp Bot 201 1, 62(8):2521-2569. 

32. Lucker J, Laszczak M, Smith D, Lund ST: Generation of a predicted protein 
database from EST data and application to iTRAQ analyses in grape {Vitis 
vinifera cv. Cabernet Sauvignon) berries at ripening initiation. 

BMC Genomics 2009, 10:1-17. 

33. Peng FY, Reid KE, Liao N, Schlosser J, Lijavetzky D, Holt R, Martinez Zapater 
JM, Jones S, Marra M, Bohlmann J, Lund ST: Generation of ESTs in Vitis 
vinifera wine grape (Cabernet Sauvignon) and table grape (Muscat 
Hamburg) and discovery of new candidate genes with potential roles in 
berry development. Gene 2007, 402(1 -2):40-50. 

34. Jaillon 0, Aury JM, Noel B, et al\ The French-Italian public consortium for 
grapevine genome characterization: the grapevine genome sequence 
suggests ancestral hexaploidization in major angiosperm phyla. 
Nature 2007, 449:463-467. 

35. Velasco R, Zharkikh A, Troggio M, et al\ A high quality draft consensus 
sequence of the genome of a heterozygous grapevine variety. PLoS ONE 
2007, 2:e1326. 

36. Martinez-Esteso MJ, Casado-Vela J, Selles-Marchart S, Elortza F, Pedreho MA, 
Bru-Martinez R: iTRAQ-based profiling of grape berry exocarp proteins 
during ripening using a parallel mass spectrometric method. Mol Biosyst 
2011,7(3):749-65. 

37. Ross PL, Huang YN, Marchese JN, Williamson B, Parker K, Hattan S, 
Khainovski N, Pillai S, Dey S, Daniels S, Purkayastha S, Juhasz P, Martin S, 
Bartlet-Jones M, He F, Jacobson A, Pappin DJ: Multiplexed protein 
quantitation in Saccharomyces cerevisiae using amine-reactive isobaric 
tagging reagents. Mol Cell Proteomics 2004, 3(1 2):1 154-69. 

38. Ong SE, Mann M: Mass spectrometry-based proteomics turns 
quantitative. Not Chem Biol 2005, 1:252-262. 

39. Casado-Vela J, Martinez-Esteso MJ, Rodriguez E, Borras E, Elortza F, Bru 
Martinez R: iTRAQ-based quantitative analysis of protein mixtures with 
large fold change and dynamic range. Proteomics 2010, 10:343-347. 

40. Pierce A, Unwin RD, Evans CA, Griffiths S, Carney L, Zhang L, Jaworska E, Lee 
CF, Blinco D, Okoniewski MJ, Miller CJ, Bitton DA, Spooncer E, Whetton AD: 
Eight-channel iTRAQ enables comparison of the activity of six 
leukemogenic tyrosine kinases. Mol Cell Proteomics 2008, 7(5):853-863. 

41 . Grimplet J, Deluc LG, Tillett RL, Wheatley MD, Schlauch KA, Cramer GR, 
Cushman JC: Tissue-specific mRNA expression profiling in grape berry 
tissues. BMC Genomics 2007, 8:187. 

42. Negri AS, Prinsi B, Rossoni M, Failla 0, Scienza A, Cocucci M, Espen L: 
Proteome changes in the skin of the grape cultivar Barbera among 
different stages of ripening. BMC Genomics 2008, 9:378. 

43. Gotz S, Garcia-Gomez JM, Terol J, Williams TM, Nueda MJ, Robles M, Talon 
M, Dopazo J, Conesa A: High-throughput functional annotation and data 
mining with the Blast2G0 suite. Nucleic Acids Res 2008, 36:3420-3435. 

44. da Silva EG, landolino A, Al-Kayal F, Bohlmann MC, Cushman MA, Lim H, 
Ergul A, Figueroa R, Kabuloglu EK, Osborne C, Rowe J, Tattersall E, Leslie A, 
Xu J, Baek JM, Cramer JR, Cushman JC, Cook DR: Characterization the 



grape transcriptome: analysis of expressed sequence tags from multiple 
Vitis species and development of a compendium of gene expression 
during berry development. Plant Physiol 2005, 239:574-597. 

45. Bluthgen N, Brand K, Cajavec B, Swat M, Herzel H, Beule D: Biological 
profiling of gene groups utilizing gene ontology. Genome Inform 2005, 
16:106-115. 

46. Storey R: Potassium localization in the grape berry pericarp by energy- 
dispersive X-ray microanalysis. Am J Enol Vitic 1987, 38:301-309. 

47. Terrier N, Francois-Xavier S, Ageorges A, Romieu C: Changes in acidity and 
in proton transport at the tonoplast of grape berries during 
development. Planta 2001, 213:20-28. 

48. Terrier N, Romieu C: Inhibition of vacuolar proton pumps by ethanol impairs 
grape berry compartmentation. AustJ Grape Wine Res 1998, 4:39-45. 

49. Shiratake K, Goto H, Maeshima M, Yamaki S: Fluctuations in the protein 
levels of H^-pumps and water channels of the tonoplast and plasma 
membrane during grape berry development. J Japan Soc Hort Sci 2001, 
70(3):287-293. 

50. Zhang X-Y, Wang X-L, Wang X-F, Xia G-H, Pan Q-H, Fan R-C, Wu F-Q, Yu X-C 
, Zhang D-P: A shift of phloem unloading from symplasmic to 
apoplasmic pathway is involved in developmental onset of ripening in 
grape berry. Plant Physiol 2006, 142:220-232. 

51. Davies JM, Poole RJ: Sanders D:the computed free energy change of 
hydrolysis of inorganic pyrophosphate and ATP: apparent significance 
for inorganic-pyrophosphate-driven reactions of intermediary 
metabolism. BBA- Bioenergetics 1993, 1 141(1 ):29-36. 

52. Plaxton WC: The organization and regulation of plant glycolysis. Annu Rev 
Plant Phys 1996, 47:185-214. 

53. Kennedy JA, Saucier C, Glories Y: Grape and wine phenolics: history and 
perspective. Am J Enol Vitic 2006, 57:239-248. 

54. Hall D, De Luca V: Mesocarp localization of a bi-functional resveratrol/ 
hydroxycinnamic acid glucosyltransferase of Concord grape {Vitis 
labrusca). Plant J 2007, 49:579-591. 

55. Kennedy JA, Matthews MA, Waterhouse AL: Changes in grape seed 
polyphenols during fruit ripening. Phytochemistry 2000, 55:77-85. 

56. Kennedy JA, Hayasaka Y, Vidal S, Waters EJ, Jones GP: Composition of 
grape skin proanthocyanidins at different stages of berry development. 
J Agr Food Chem 2001 , 49:5348-5355. 

57. Downey MO, Harvey JS, Robinson SP: Analysis of tannins in seeds and 
skins of Shiraz grapes throughout berry development. Aust J Grape Wine 
Res 2003, 9:15-27. 

58. Lima MRM, Ferreres F, Dias ACP: Response of Vitis vinifera cell cultures to 
Phaeomoniella chlamydospora: changes in phenolic production, 
oxidative state and expression of defence-related genes. Eur J Plant 

Pathol 2012, 132(1):133-146. 

59. Bogs J, Downey MO, Harvey JS, Ashton AR, Tanner GJ, Robinson SP: 
Proanthocyanidin synthesis and expression of genes encoding 
leucoanthocyanidin reductase in developing grape berries and 
grapevine leaves. Plant Physiol 2005, 139:652-663. 

60. Zamboni A, Di CarIi M, Guzzo F, Stocchero M, Zenoni S, Ferrarini A, Tononi 
P, Toffali K, Desiderio A, Lilley KS, Pe ME, Benvenuto E, Delledonne M, 
Pezzotti M: Identification of putative stage-specific grapevine berry 
biomarkers and omics data integration into networks. Plant Physiol 2010, 
154:1439-1459. 

61. Marrs KA: The functions and regulation of glutathione S-transferases in 

plants. Annu Rev Plant Phys 1996, 47:127-158. 

62. Mueller LA, Walbot V: Models for vacuolar sequestration of anthocyanins. In 
Regulation of phytochemicals by molecular techniques. 35th edition. Edited by 
Romeo JT Saunders JA, Matthews BE. New York: Pergamon; 2001:297-312. 

63. Terrier N, Olle D, Verries C, Cheynier V: Biochemical and molecular aspects 
of flavan-3-ol synthesis during berry development. In Grapevine molecular 
physiology and biotechnology. Edited by Roubelakis-Angelakis KA. New York: 
Springer; 2009:365-388. 

64. Ageorges A, Fernandez L, Vialet S, Merdinoglu D, Terrier N, Romieu C: Four 
specific isogenes of the anthocyanin metabolic pathway are 
systematically co-expressed with the red colour of grape berries. 

Plant Sci 2006, 170:372-383. 

65. Martinez-Esteso MJ, Selles-Marchart S, Vera-Urbina JC, Pedreho MA, 
Bru-Martinez R: DIGE analysis of proteome changes accompanying large 
resveratrol production by grapevine {Vitis vinifera cv. Gamay) cell 
cultures in response to methyl-(3-cyclodextrin and methyl jasmonate 
elicitors. J Proteomics 201 1, 74(8):1 421 -1436. 



Martinez-Esteso et al. BMC Plant Biology 201 3, 1 3:1 67 Page 20 of 20 

http://www.bionnedcentral.conn/1471-2229/13/167 



66. 



67. 



69. 
70. 
71. 

72. 
73. 

74. 

75. 

76. 
77. 



79. 



81. 



85. 



87. 



Bilang J, Macdonald H, King PJ, Stum A: A soluble auxin-binding protein 
from Hyoscyamus muticus is a glutathione S-transferase. Plant Physiol 

1993, 102:29-34. 

Gonneau J, Mornet R, Laloue M: A Nicotiniana plumbaginifolia protein 
labeled with na azido cytokinin agonist is a glutathione S-transferase. 

Physiol Plant 1998, 103:1 14-124. 

Koeduka T, Louie GV, Oriova I, Kish CM, Ibdah M, Will<erson CG, Bowman 
ME, Baiga TJ, Noel JP, Dudareva N, Pichersl<y E: The multiple 
phenylpropene synthases in both Clarkia breweri and Petunia hybrida 
represent two distinct protein lineages. Plant J 2008, 54(3):362-374. 
Fenoll J, Manso A, Hellin P, Ruiz L, Fiores P: Changes in the aromatic 
composition of the Vitis vinifera grape Muscat Hamburg during ripening. 
Food Chem 2009, 1 14:420-428. 

Pandey RM, Farmahan HL: Changes in the rate of photosynthesis and 

respiration in leaves and berries of Vitis vinifera grapevines at various 

stages of berry development. Vitis 1977, 16:106-1 11. 

Famiani F, Walker RP, Tecsi L, Chen Z-H, Proietti P, Leegood RC: An 

immunohistochemical study of the compartmentation of metabolism 

during the development of grape {Vitis vinifera L.) berries. J Exp Bot 2000, 

51:675-683. 

Ollat N, Gaudillere J: Carbon balance in developing grapevine berries. 
Acta Hon. (ISHS) 2000, 526:345-350. 
Grace S, Pace R, Wydrzynski T: Formation and decay of 
monodehydroascorbate radicals in illuminated thylakoids as determined 
by EPR spectroscopy. Biochinn Biophys Acta 1995, 1229:155-165. 
Yang EJ, Oh YA, Lee ES, Park AR, Cho SK, Yoo YJ, Park OK: Oxygen-evolving 
enhancer protein 2 is phosphorylated by glycine-rich protein 3/wall- 
associated kinase 1 in Arabidopsis. Biochem Bioph Res Co 2003, 305:862-868. 
Yang Y, Sulpice R, Himmelbach A, Minhard M, Christmann A, Grill E: Fibrillin 
expression is regulated by abscisic acid response regulators and is 
involved in abscisic acid-mediated photoprotection. Proc Natl Acad Sci 
USA 2006, 203(1 5):6061 -6066. 

Campanella M, Parker N, Tan CH, Hall AM, Duchen MR: IF1: setting the 
place of the F1 FO-ATP synthase. Trends biochem Sci 2009, 34(7):343-350. 
Sweetman C, Deluc LG, Cramer GR, Ford CM, Soole KL: Regulation of 
malate metabolism in grape berry and other developing fruits. 

Phytochennistry 2009, 70:1 329-1 344. 

Pilati S, Perazzolli M, Malossini A, Cestaro A, Dematte L, Fontana P, Dal Ri A, 
Viola R, Velasco R, Moser C: Genome-wide transcriptional analysis of 
grapevine berry ripening reveals a set of genes similarly modulated 
during three seasons and the occurrence of an oxidative burst at 
veraison. BMC Genomics 2007, 8:428. 

Tesniere C, Verries C: Molecular cloning and expression of c-DNAs 
encoding alcohol dehydrogenases from Vitis vinifera L. during berry 
development. Plant Sci 2000, 157:77-88. 

Ponce-Valadez M, Watkins CB: Fermentation and malate metabolism in 
response to elevated C02 concentrations in two strawberry cultivars. 

Physiol Plant 2008, 1 34:1 21 -1 33. 

ChareonthiphakornN WD: Golan-GoldhirshA Wititsuwannakul R: 
purification and characterization of NAD(P)H quinone reductase from 
the latex of Hevea brasiliensis Mull.-Arg. (Euphorbiaceae). 

Phytochemistry 2002, 61 (2):1 23-1 28. 

Deller S, Macheroux P, Sollner S: Flavin-dependent quinine reductases. 

Cell Mol Life Sci 2008, 65:1 41 -1 60. 

Sollner S, Schober M, Wagner A, Prem A, Lorkova L, Palfey BA, Grol M, 
Macheroux P: Quinone reductase acts as a redox switch of the 20S yeast 
proteasome. EMBO J 2009, 10:65-70. 

Agasse A, Vignault C, Kappel C, Conde C, Geros H, Delrot S: Sugar transport and 
sugar sensing in grape. In Grapevine molecular physiology and biotechnology 
Edited by Roubelakis-Angelakis KA. New York: Springer; 2009:105-139. 
Dali N, Michaud D, Yelle S: Evidence for the involvement of sucrose 
phosphate synthase in the pathway of sugar accumulation in sucrose- 
accumulating tomato fruits. Plant Physiol 1992, 99:434-438. 
D'Aoust MA, Yelle S, Nguyen-Quoc B: Antisense inhibition of tomato fruit 
sucrose synthase decreases fruit setting and the sucrose unloading 
capacity of young fruit. Plant Cell 1 999, 1 1 :2407-241 8. 
Nguyen-Quoc B, Foyer CH: A role of 'futile cycles' involving invertase and 
sucrose synthase in sucrose metabolism of tomato fruit. J Exp Bot 2001, 
52:881-889. 



90. 

91. 
92. 

93. 

94. 

95. 

96. 
97. 



99. 
100. 

101. 

102. 

103. 
104. 



Takayanagi T, Yokotsuka K: Relationship between sucrose accumulation 
and sucrose-metabolizing enzymes in developing grapes. Am J Enol Vitic 

1997, 48:403-407. 

Koch KE, Nolte KD, Duke ER, McCarty DR, Avigne WT: Sugar levels 
modulate differential expression of maize sucrose synthase genes. 

Plant Cein992, 4:59-69. 

Steffan H, Rapp A, Ullemeyer H, Kupfer G: Uber den reifeabhangigen 
Saure-Zuker-Stoffwechsel bei beeren von V. Vinifera-Sorten, untersucht 
mit ^"^C Verbindungen. Vitis 1975, 14:181-189. 

Coombe BG: Distribution of solutes within the developing grape berry in 
relation to its morphology. Am J Enol Vitic 1987, 38:120-127. 
Or E, Baybik J, Sadka A, Saks Y: Isolation of mitochondrial malate 
dehydrogenase and phosphoenolpyruvate carboxylase cDNA clones 
from grape berries and analysis of their expression pattern throughout 
berry development. J Plant Physiol 2000, 157:527-534. 
Hawker JS: Changes in the activities of malic enzyme, malate 
dehydrogenase, phosphopyruvate carboxylase and pyruvate 
decarboxylase during the development of a non-climacteric fruit (the 
grape). Phytochemistry 1969, 8:19-23. 
Diakou P, Svanella L, Raymond P, Gaudillere JP, Moing A: 
Phosphoenolpyruvate carboxylase during grape berry development: 
protein level, enzyme activity and regulation. Aust J Plant Physiol 2000, 
27:221-229. 

Ruffner HP, Hawker JS, Hale CR: Temperature and enzymic control of 
malate metabolism in berries of Vitis vinifera. Phytochemistry 1 976, 

15:1877-1880. 

Ruffner HP, Kliewer WM: Phosphoenolpyruvate carboxykinase activity in 
grape berries. Plant Physiol 1 975, 56:67-71 . 

Taureilles-Saurel C, Romieu CG, Robin JP, Flanzy C: Grape (Vitis vinifera L.) 
malate dehydrogenase: II: Characterization of the major mitochondrial and 
cytosolic isoforms and their role in ripening. Am J Enol Vitic 1 995, 46:29-36. 
Ruffner HP, Possner D, Brem S, Past DM: The physiological role of malic 
enzyme in grape ripening. Planta 1984, 160:444-448. 
Davies DD: The fine control of cytosolic pH. Physiol Plant 1986, 67:702-706. 
Ruffner HP: Metabolism of tartaric and malic acids in Vitis: a review - Part 
B. Vitis 1982, 21:346-358. 

Peruffo ADB, Pallavicini C: Enzymatic changes associated with ripening of 

grape berries. J Sci Food Agric 1 975, 26:559-566. 

Bradford MM: A rapid and sensitive method for the quantitation of 

microgram quantities of protein utilizing the principle of protein-dye 

binding. Anal Biochem 1976, 72:248-254. 

Stum A, Quackenbush J, Trajanoski Z: Genesis: cluster analysis of 

microarray data. Bioinformatics 2002, 18(1):207-208. 

Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam 

H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, 

Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics 2007, 

23:2947-2948. 



doi:1 0.1 1 86/1 471 -2229-1 3-1 67 

Cite this article as: Martinez-Esteso et al.: iTRAQ-based protein profiling 
provides insights into the central metabolism changes driving grape 
berry development and ripening. BMC Plant Biology 2013 13:167. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



